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Abstract 

We determine all spacetime flows which lead to a Hamiltonian dynam¬ 
ics in the space of general relativistic initial data sets with asymptoti¬ 
cally Kerr-de Sitter ends. The corresponding Hamiltonians are calculated. 
Some implications for black-hole thermodynamics are pointed out. 
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1 Introduction 

There is growing astrophysical evidence that spacetimes with positive cosmo¬ 
logical constant provide physically correct models for cosmology. Large fam¬ 
ilies of such non-compact, vacuum, general-relativistic models have been con¬ 
structed (see [15,16,19, 20, 22, 23,41] and references therein). In particular we 
have now large classes of initial data sets with one or more ends of cylindrical 
type, in which the metric becomes periodic when one recedes to infinity along 
half-cylinders. This then raises the question, how to define the total mass, or 
energy, of such configurations. We have answered this question for asymptot¬ 
ically Schwarzschild-de Sitter metrics in [12] using a Hamiltonian formalism, 
and the aim of this paper is to generalize the analysis there to metrics with 
asymptotically Kerr-de Sitter ends. 

Consider, thus, the variational identity associated with the motion of a 
hypersurface in a vacuum space-time in the direction of a vector field X, Equa¬ 
tion (2.3) below. The requirement of convergence of the volume integral there, 
when the boundary recedes to infinity on an asymptotically periodic end, leads 
to the requirement that the Lie derivative of the metric in the direction of X 
approaches zero asymptotically, at least in an integral sense. Barring the case 
where X is everywhere spacelike at large distances, one expects that this will 
only happen for metrics which are asymptotic to Kerr-de Sitter metrics. Hence 
our interest in asymptotically Kerr-de Sitter metrics. We note that a large fam¬ 
ily of non-trivial vacuum spacetimes which are exactly Kerr-de Sitter outside 
the domain of dependence of a compact set has been constructed in [23]. 
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Now, in the coordinate system of (3.1) below, all Killing vectors of Kerr-de 
Sitter metrics with rotation parameter a 7 ^ 0 take the form 

X = xdt + f3d^, ( 1 . 1 ) 

where x and /3 do not depend upon the spacetime coordinates, but might depend 
upon the solution at hand. 

The vector field dip is singled-out by the requirement that all its orbits are 
27r-periodic. Not unexpectedly, we check that the spacetime flow of dp generates 
a Hamiltonian flow on the space of metrics. The corresponding Hamiltonian is 
usually interpreted as the total angular momentum of the solution. 

In the case where the rotation parameter a of the metric vanishes, the vector 
dt is singled-out as the unique, up to a multiplicative constant. Killing vector 
field orthogonal everywhere to dp. We showed in [12] that the spacetime flow of 
dt generates a Hamiltonian dynamical system in the space of gravitational ini¬ 
tial data with asymptotically Schwarzschild-de Sitter ends, with the associated 
Hamiltonian equal to the mass parameter m. 

When a 7 ^ 0, there does not seem an obvious choice for a second preferred 
Killing vector. A natural question arises, which of the Killing vector fields of 
the form (1.1) generate a Hamiltonian dynamical system on the space of asymp¬ 
totically Kerr-de Sitter vacuum metrics. We give an exhaustive description of 
such vector fields in Theorem 5.1 below. This is the main result of this work. 
The result requires a detailed analysis of the variational identities arising in 
this context, to be found in Section 2 which occupies a significant part of this 
paper. We believe this analysis has interest on its own. 

Much to our surprise, the vector field dt of the coordinates of (3.1) is 
not associated with a Hamiltonian flow. However, we show that the vector 
field Y^l -|- a?K/2,dt is, the corresponding Hamiltonian being equal to m(l -|- 
a^A/3)“^/^. Some further explicit examples of interest of Hamiltonian flows 
are given. 

As such, the simplest asymptotically Kerr-de Sitter spacetime is the Kerr-de 
Sitter spacetime itself. The variational identities that result from a Hamilto¬ 
nian analysis are then known in the literature under the name of “black hole 
thermodynamics”. The ambiguity in the choice of the energy, equivalently, in 
the choice of the Hamiltonian vector field, leads to ambiguities in the definition 
of thermodynamical variables. Here we point out that in the Kerr-de Sitter 
spacetimes one can define preferred Killing vector fields X by requiring that X 
is tangent to some bifurcate Killing horizon, and has surface gravity equal to 
one there. (This enables one to define e.g. the rotation velocities of the remain¬ 
ing horizons with respect to the selected one.) It turns out that these vector 
fields are Hamiltonian, with Hamiltonian equal the area of the selected horizon. 
The choice of “total energy” as the value of the associated Hamiltonian, leads 
to a set of geometrically unique thermodynamical identities. 

We note the recent analysis of [30], where thermodynamical instabilities are 
tied to instabilities of associated black string solutions. We believe that our 
analysis provides a good starting point for similar considerations for solutions 
with a positive cosmological constant, and hope to be able to address this issue 
in a near future. 
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We refer the reader to [1, 4-6, 39] and references therein for alternative ap¬ 
proaches to a definition of mass in the presence of a positive cosmological con¬ 
stant. In Appendix F we shortly discuss a definition using conformal Killing- 
Yano tensors. We note that in space-times with a non-zero cosmological con¬ 
stant the Komar integral does not provide a surface-independent integrand ex¬ 
cept for the angular-momentum of hypersurfaces asymptotically invariant under 
rotations. 

2 The variational formula 

We wish to define a Hamiltonian dynamical system on a set of Lorentzian 
metrics on an (3 -|- l)-dimensional manifold ^^ assuming the existence of a 
spacelike hypersurface 5^ C ^ on which the metric asymptotes to a Kerr-de 
Sitter metric as one recedes to infinity along an end of . This is a special 
case of a more general construction which can be done in all space-dimensions 
n > 2 and proceeds as follows: Given a Lorentzian metric g on ^ we choose 
an n-dimensional manifold S and a one-parameter family of embeddings 

M X S 9 (t, x) ip{t, x) G ^ . (2.1) 

Note that while for all t G M the maps Y 9 x i— >• '0(t, x) G ^ are assumed to be 
embeddings of Y, we do not assume that the whole map (2.1) is an embedding. 
Thus, the hypersurfaces 


■= {V'(Lx) , X G Y} 

do not need to form a foliation, they are allowed to cross each other, etc. 

We emphasise that the maps V’ can depend upon the metric g under con¬ 
sideration. Indeed, this will be the case in our main application here, namely 
the description of a family of Hamiltonian dynamical systems in the space of 
asymptotically Kerr-de Sitter metrics. 

To continue, we choose a compact volume K C Y with a smooth boundary 
dV. The equations will be analysed on V before taking an exhaustion of Y and 
passing to the limit. We set 

r = V'(0,K). 

Along each hypersurface we define a spacetime vector field X by the 
formula 

Y := . (2.2) 

Note that this defines a vector field on a neighbourhood of =5^ ;= if if) is 
a diffeomorphism near {0} x Y but, in general, we will only have a spacetime 
vector field defined along for each t. 

2.1 Adapted coordinates 

In [33] it was assumed that X is actually a vector field on ^, which is moreover 
supposed to be timelike and non-vanishing. Then, the spacetime domain 

VL-.= ut{G^{y)} 
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was considered, where by {Qf} we denote the (local) group of diffeomorphisms 
generated by the field X. The boundary defined as dQ := VJt{G^{dV)} was, 
therefore, a smooth timelike submanifold, with 'V := 5^. Let (x^), k = 

1,... ,n, be a coordinate system on S. In [33] local coordinates on ^ were 
chosen so that 

= t}, and y*(V’(t,x)) = x*(x), 

whence X = d^. Using this particular coordinate system, the following vari¬ 
ational formula has been proved in [33] in dimension 3-1-1, for a Lorentzian 
metric interacting with a matter field tp: 

0 = ^ [ (P^'Sgki - + [ {pSp - pdp) + - [ {XSa - aSX) 

Jv^ ^ Jv 7 Jdv 

+ ^ [ {2iy5Q - 2iy^5ClA + vQ^^^gAB) + [ pSp (2.3) 

^7 JdV JdV 

with 7 = Stt when n = 3 (see, e.g., [12, Appendix D] for other dimensions). 
Here, gki and P^^ describe gravitational Cauchy data on V, i.e. the pull-back to 
V of the metric induced by the space-time metric on P and of the ADM mo¬ 
mentum, respectively, whereas p and p represent symbolically the corresponding 
Cauchy data for the matter fields, if any. “Dot” stands for the time derivative 
dt. As explained in detail below or in [33], the Q’s describe various components 
of the extrinsic curvature of the world tube dGl, whereas u, and gAB en¬ 
code the metric induced on dXl using an ((n — 1) -|- 1) decomposition. Finally, 
A = Vdet gAB is the volume form on dV whereas a denotes the hyperbolic angle 
between XP and dGl. For a Hamiltonian flow the second line of (2.3) equals the 
variation of the Hamiltonian, which is at the origin of the occurrence of the 7 
factors in the formula. 

As such, in [28] the formula was generalized to the case when X is spacelike. 
Next, in [24] the held was allowed to be null-like. Finally, it has been pointed 
out in [12] that the formula remains true for vacuum metrics, possibly with a 
cosmological constant, in any space-dimension n > 2, with a constant 7 which 
depends upon dimension. 

Our hrst aim is to remove the coordinate condition X = do, and the hypoth¬ 
esis that A is a vector held on ./#, and allow for non vanishing variation 6X^ of 
the held X. The latter is motivated by the fact, that imposing coordinate con¬ 
ditions leads typically to vector helds X which depend upon the conhguration 
of the held, so varying the latter necessarily requires varying the former. This 
issue will indeed have to be addressed in our analysis below of asymptotically 
Kerr-de Sitter spacetimes. 

Now, on each V) := {t} x U we have a non-degenerate n-dimensional metric g 
and the ADM momentum P, dehned via the V'-pull-back of the corresponding 
objects from 'ip{{t} x U) G Given a coordinate system (x^) on V, we 

denote components of g and P, together with their derivatives with respect to 
the parameter t, as gki, P^\ gu and respectively. On each boundary dVt 
we have the induced (n — l)-metric, which we describe by its components gAB- 
For computational purposes it is useful to choose the last coordinate x” in such 
a way that it is constant on dV, and then the collection (x"^), A = 1,..., n — 1, 
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provides a coordinate chart on the boundary. By A := \J det gAB we denote the 
(n — l)-volume density on the boundary. 

At each point of d'f we decompose the field X into the part tangent to 91^, 
which we denote by and the part A-*-, orthogonal to dX. We have, 

therefore, A = A-*- + A^ and we define 

^ '■= 

where the + sign is taken if A is timelike and ” if A is spacelike. At points 
at which A-*- is non-zero and not null we define the unit vector N := ^A-*-, so 
that there we have; 

X = uN + v^Oa . (2.4) 

We will use the same symbols and v for the corresponding pull-backs to dV. 

To define the remaining objects appearing in (2.3), on dX consider the 
two-dimensional tangent plane in orthogonal to TdY (the plane may be 
identified with the two-dimensional Minkowski space) and use the following 
four normalized vectors: N, M - orthogonal to N, m - tangent to di¬ 
rected outwards and, finally, n - orthogonal to m, directed in the future. See 
Figure 2.1. 

By “a” we denote the “hyperbolic angle” between N and n, defined as 
follows: 

J arsinh(N | m) for A-*- timelike, , . 

\ sgn(N I m)arcosh(N | m) for A"*- spacelike. 


Of course, the definition is meaningful only for non-vanishing which is the 
case if we assume that A-*- is everywhere timelike, or everywhere spacelike non¬ 
vanishing everywhere. The hypothesis that v has no zeros is needed for the 
derivation of the elegant formula (2.3). However, we stress that this hypothesis 
will not be needed in our main results below. 

In order to define the remaining objects Q, Qa and we consider the 
hypersurface cAl = {Q^{dV)} C parameterized by coordinates (x“) = 
= {t,x^), and the ADM version of its extrinsic curvature: 

:= Vldet^cdl (Tr'-T“'), (2.6) 

= (2-7) 

where is the n-dimensional inverse with respect to the induced metric gab 
on with L = L°'^gah- 
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Assuming again that v has no zeros, we set: 


Q = 

Qa = Q\ = Q^^gaA. 

qab ^ , 

V V ^ ' 

where we use (n — l)-dimensional inverse of the metric qab on dV. 

We can also reformulate these definitions in terms of the extrinsic curvature 
of d'Y^ namely; the torsion covector 

£a:=(VaN|M), (2.8) 

and the extrinsic curvature tensor in direction of M 

kAB = kABiM) := (Va^bIM) , (2.9) 

and its trace 

'^AB 

k = g kAB ■ (2.10) 

Equivalently, for any pair (y, Z) of vector fields tangent to dY it holds that 
k{Y,Z) = (VyZ|M). We have 

Q = Xk, (2.11) 

Qa = XiA, (2.12) 

= X(^k^^-g^^k + g^^s^ , (2.13) 

where by s we denote the following “acceleration scalar” 

s = (V^NIM) =-(M| V^M). (2.14) 

This completes the list of geometric objects used in (2.3). 

In Appendix C we prove the following;^ 

Theorem 2.1 Formula (2.3) is valid for V which are either spacelike every¬ 
where or timelike everywhere, for any veetor field X defined along V such that 
X is nowhere vanishing and everywhere transverse to the image X of V in 
spacetime. The vector field X is allowed to depend upon the metric and its 
derivatives (hence with non-vanishing variation 5X>^) provided that X is kept 
fixed in spacetime. 


^More precisely, in Appendix C we prove the second part of our theorem, namely that the 
vector field X is allowed to depend upon the field configuration without changing (2.3). The 
first part of the theorem is implicit in the considerations in [33], and follows explicitly from 
the remaining analysis in the current work in any case. 
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2.2 General case 


As made clear in Theorem 2.1, the analysis of (2.3) so far assumed that the 
field A is a) a spacetime vector field b) which is metric-independent and c) 
is transverse to 'f' with spacelike or timelike. However, in our analysis of 
asymptotically Kerr-de Sitter metrics we need to allow vector fields X which 
depend upon the metric, which might vanish and/or change type, and in fact 
hypersurfaces Y which might be either spacelike everywhere or timelike every¬ 
where. We will remove all these unwanted conditions. This forces us to revisit 
the calculations leading to (2.3) for general vectors X. 

The essential difference between our final formula below and that in [33] is, 
that here the geometric quantities appearing on dV are not necessarily those 
associated with a hypersurface obtained by flowing dX along X, as this set will 
not form a hypersurface in general, but those associated with a hypersurface 
{y” = 0} obtained by choosing some coordinate system near X in which X = 
{y^ = 0}, dX C {y” = 0}. When X is everywhere transverse to X, we can 
choose coordinates so that X = d/dy^, in which case all objects appearing in 
this section coincide with those of [33], as described in Section 2.1. However, 
transversality will not be assumed in this section. 

We emphasise that even in the transverse case we do not need to make 
the choice X = d/dy^ (or X proportional to d/dy^, which leads to the same 
geometric objects). There is indeed a freedom here, which is equivalent to the 
choice of the vector field d/dy^ at the boundary dY. Every such choice will 
lead to the variational formula (2.35) below, with different geometric quantities 
appearing there at dY for two not-everywhere-colinear-at-9A choices of d/dy^. 
A natural possibility is to choose d/dy^ to be normal to Y, but other choices 
might be more convenient depending upon the problem at hand. 

In any case, we stress once again that the index 0 does not indicate the 
variable t associated with Hamiltonian flow, but an auxiliary coordinate y®. 

To continue, consider the Lagrangian for vacuum Einstein’s equations with 
a cosmological constant A; 

~ > (2-15) 

where is expressed in terms of a metric and a symmetric connec¬ 

tion T^^. The multiplicative coefficient I/IGtt is physically motivated only in 
dimension 3-1-1, using units where G = 1. Following [33] we define 

Consider a one-parameter family of field configurations A (y^,y(A), r^j^(A)) 
and define a “variation” S as the operation 


5 ;= 


d 

dX 


A=0 


thus 5g = ^|a= 0 ) etc. The following identity can be easily verified: 
6L = + 



=:S’P-^ --^A 

= - (V,P/"") + d^ (p/"".5r;:,) , (2.16) 

where 

p^^^un ^ ^ ( 2 . 17 ) 

Due to the identity (2.16), the (vacuum) Euler-Lagrange equations, 


0 = 


5L 



are equivalent to 

SL = d^ (p/"".5r;:. 

To simplify notation we set 




(2.18) 

(2.19) 


( 2 . 20 ) 

which enables us to rewrite identity (2.16) in yet another, equivalent form. 

SL - = dx [n^’^SAf^,) . (2.21) 

This equation holds regardless of whether or not the metric g satisfies itself the 
vacuum field equations or the connection fulfills the metricity condition. 

Now, each field configuration A i—)■ g{X) comes with its own family of maps 

(A,t) 1-^ 'ipxit,-) 

as at the beginning of Section 2, where each •) : S i—)■ ^ is an embedding 
of S D E into the associated spacetime Let 


A I—)■ X(A) :— 


be the one-parameter family of spacetime vector fields defined along ipxiit} x S) 
associated with ipx- At given t, the field X(A) can be extended to a smooth 
vector field defined in a spacetime neighborhood of 'i/’A({^} P) in many ways. 
If the V'a’s are local diffeormorphisms near {t} x S, then ip\dt defines directly 
such a vector field, and we will make this choice. Otherwise we choose any 
smooth extension; one of the important outcomes of our calculations below is 
to show that the result does not depend upon this choice. 

To derive the “Hamilton-type” variational identity, (2.27) below, we rewrite 
the Lie derivative of the connection, 

, ( 2 . 22 ) 


in terms of A (compare [33, 34]) 


{lU 


{lU 


VW\ [ 

Ibvr \ 
1 


h 


1 9 a 

167ri ^ 


V/,(v'^x“ - v"a:^) + 2P“^a:‘"| 
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. (2.23) 



This leads to the identity 


+d^ (2.24) 

+L5X^ + 7r^"[jSfx, , 

keeping in mind that we allow X to depend upon the field configuration, and 
where we do not assume that the vacuum Einstein equations are satisfied. How¬ 
ever, we have assumed that the connection T is metric; equivalently, the field 
defined in (2.16) vanishes. Using (2.23), we find 

- X^L + 

[v 1^(V^X^ - V^X^)] + ( 2v1^i?\X^" - 167rX^L )| . (2.25) 

=-327rX'^£’^,^ 


Recall that 


The identity \^x , <5] 
follows: 


5X 


_ ^ 
“ lx 


—A£^x enables one to rewrite the last line in (2.24) as 


UX^ - = 

[v1^(V^<5X^ - V^JX^)] - 327r£:\,5X'^| . (2.26) 


We conclude that for solutions of the field equations the term — 

^xA^j^dx'^'' in (2.24) is the divergence of a bivector density: 

- {^xTrnSA^^u = -2X^6£\ + 

+Y^5m{<5[v1^(v^x^ - v^x^)]| + d, [x\^^^5Ai, - 

\^M{^^^X^ - V^<5X^)1 . (2.27) 

Let (y°,y*) be local coordinates near X = V'o({0} x U) = V'dO} x V) so 
that tp is constant on X. The y^’s here will always be local coordinates on ^, 
with = djdy^^ while {t, x*) will always denote local coordinates on M x S. In 
particular we emphasise that do = d/dy^ will not be equal to d/dt in general. 

From now on we assume that X is spacelike everywhere or timelike ev¬ 
erywhere; equivalently, is nowhere vanishing on X. As such, studies of 
Hamiltonian dynamics assume that X is spacelike. However, for the analysis 
in Section 4 the variational identities that we are about to derive will also be 
needed for X’s which are timelike. 
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In coordinates such that Y C {y^ = const}, the integral of the left-hand 
side of (2.27) over the n-dimensional manifold 'f' reads 


ly 


- i^xTr^SA' 




dSA 


Jy 


(2.28) 


We wish to rewrite this formula in terms of the usual initial data on For 
this, set 




dy^ 


A 

dyk 


where Y := X^-—r is tangent to Y. 


Let / be some field on spacetime. We shall denote by ^xf the usual 
Lie-derivative operator on a manifold, and by £x/ the restriction of ^xf to 
Y. Note that jSfyF coincides with CyF for vector fields Y tangent to Y and 
geometric fields F on F, and both notations will be used interchangeably in 
such cases. See Appendix B for some explicit expressions. 

We have the following: 


Theorem 2.2 


1. Suppose that X = X^Oq on F, then: 

[ (£x9kiSP’^'-FxP^^Sgki)+2[ {CxaSX - CxXha) 
Jy ^ ' Jdy 

- [ A°(2z/5Q - 2iy^6QA + nQ^^hgAB) 

Jay 


where 


= IGtt / aO {-25S\ + , 

Jy 


Cx09,gki := X^dogki + NkdiX^ + NidkX\ 

+ {V9 - 2f^D^N^ 

+ pkl^r _ pkr^l _ ^ 

:= A^clgA -|- Xb^OaX^ , 


(2.29) 

(2.30) 


■— X^doa— 


N~g 


,nk 


VW 


-MX^ . 


(2.31) 

(2.32) 

(2.33) 


2. Suppose that X^ = 0 X = X^di =: Y on F, then: 
(CvgkiSP^^ - FyP^^^ 


2Y^6P^k - Y^P’^^Sgki) - 327r / Y^SS^i , (2.34) 

lay ' Jy 

where Cy = Afy is the usual Lie-derivative operator on F. 
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Remarks 2.3 1. Adding (2.29) and (2.34), one obtains the variational for¬ 
mula for a general vector field X = + Y : 


ty 


(yC.x9ki^P^^ - - Px^do^^Oi) 


- / X'^{2v5q, - 2 u^6Qa + lyQ^^SgAB) 

Jar 


'dr 


2Y’^6P^k - Y^P^^ 


= IGtt / - 2X>^6rf,, (2.35) 

Jy 


with Cx for gij and obtained by adding Cxoq^ and Cy- 

2. For solutions of the field equations, (2.29) with X^ = 1 clearly coincides 
with Kijowski’s formula (2.3). 

3. We show in Appendix D below that (2.3) with the additional condition 
that Y is transversal to dP (equivalently, v has no zeros) coincides with 
(2.34). 


Proof: Since the proofs are quite lengthy and computationally intensive, we 
start with roadmaps. 

The proof of formula (2.34) proceeds as follows: 

• We integrate (2.27) over P and show that the left-hand side gives (2.71). 

• Using (2.50), (2.52) and (2.69), the boundary terms on the right-hand 
side of (2.68) lead to the formula (2.72), which is equivalent to (2.34). 

For (2.29) more work is needed: 

• We integrate (2.27) over P and show that the left-hand side splits into 
several terms as in (2.73). The terms proportional to X^ are given by 
(2.82), while the next term is a straightforward divergence. The terms 
involving d^X^ are the most problematic ones, we return to them in the 
last three steps of the calculation. 

• The divergence terms from the right-hand side of (2.27) are collected as 
a boundary integral in (2.89). The second line of this formula simplihes 
to (2.90). 

• An intermediate result is provided by (2.91), where all volume terms con¬ 
taining dkX^ are collected in (2.92). The proof is then complete when X^ 
is space-independent, as all dkX^ terms vanish under this hypothesis. 

• To analyze (2.92) one introduces an ADM-type parameterization of the 
fields. A lengthy calculation gives (2.98). 

• Using the Lie derivatives worked-out in Appendix B, one is led to (2.104). 
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• All boundary terms containing cancel out and we obtain our final 

formula (2.29). 

Let us pass now to the details of the argument. We start by recalling the 
transport formula for Christoffel symbols. Letting L denote the connection co- 
effiences in a coordinate system y, and L those transported by a diffeomorphism 
tp, it holds that 


- dy°‘ dipP dy^ 

Pi dip'^ dyl^ dy'^ dip'^ dy'^dy^ 


(2.36) 


where we have denoted by dy/dip the derivatives of the map inverse to ip. Let 
us denote by those variations of which arise from a family of maps 

'0(A) such that 0(0) is the identity. It follows from (2.36) that 






^xA% 


= 5XPd,T%-T-f,^- 

d'^SX^ 

^ dy'^dy!^ ’ 


5XPd,A-,^-A-,^ — 


d5X<^ dSX^ 85X° 

dya ^ QyP ^ Qy^ 


d‘^5X^ d^5X^ 

~^dy~fdy^ dyAdy^ ’ 

d'^X^ d‘^X^ 

dy'^ dyl^ dyA dy^ 


, dx- dx- 

-7 QyP ^ P- Qy^ 


(2.37) 


(2.38) 


(2.39) 


Now, it should be kept in mind that in all formulae above part of the 
variations of the field arise from variations of the map 0, which is allowed to 
vary. To make this clear, let us denote hy 5f those variations of the fields for 
which 0 is the identity (here the index / stands for “fields”). We thus have 


5 — Sf + 5^ . 


(2.40) 


The distinction between 5 and 5f is somewhat arbitrary for tensor fields defined 
in spacetime in a fixed coordinate system, since 5^ can always be absorbed in 
a redefinition of 6f and vice-versa. The question arises, whether this remains 
true after pull-backs to the image of 0o(0 •) have been taken. In order to clarify 
this, let us denote by gij the metrics induced by g^y{\) on the images of 0 a (L •) 


gij{X,t,x) = c/^,,(A,0a(0T)) 


dTpl{t,x) 


Then, in coordinates such that y'^(0o(O, x^)) = 0, y*(0o(O,x^)) = xL 


^f9ij (0) ■ 


gfAu{X,ipx{o,x)) 


<90(((O,x) 90^(0, x) 
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dgij{\,x) 


dX 


(2.41) 


A=0 


dxpgij{0,x') — 


d f rn , rn ^V’a(0,x) 

5a - 


dx^ 


A=0 


= ^5x9ij (2.42) 

(space-components of a spacetime Lie derivative), 


5gij{0,x) = 


d f ^aa^V’a(0>®)5V'a(0,x) 

5a - 


(9x* 


A=0 


= <J/ffij(0,x)-f (5^5ij(0,x). 


(2.43) 


Since our variations are arbitrary, for any d^gij we can redefine dfgij so that 5gij 
takes any desired value, and vice-versa. This remains true even if the variations 
are constrained to satisfy the linearised field equations (which we do not assume 
in most of our calculations, but which might be convenient for some purposes), 
since <5^-variations do indeed satisfy the linearised field equations, respectively 
the linearized constraint equations, when the fields being varied satisfy the full 
equations, respectively the constraint equations. 

From now on we will not make a distinction between gij and gij. 

A similar argument applies to Henceforth, from now on we will not 
make a distinction between 5f and 5 unless a significant ambiguity arises. 

We continue with a preliminary result which deserves highlighting; 

Proposition 2.4 The integral (2.28) depends upon the extension of X = 
off X only through boundary terms arising at dX. 


Proof: We will need the constraints implied by the identities = 0 and 

= 0: Indeed, expressing the left-hand sides in terms of and we 
obtain the following constraints; 

■^00 = -fa (4»“ + . (2.44) 

= -^(3'=’f“ + 24l>r“) ■ (2.45) 

To continue, we insert (2.39) into (2.28). We start by noting that all second- 
order derivatives there with n = n = 0 cancel out. Next, all terms 

involving doX^ and d^doX^ can be collected into 


(1) := / 4(9oXO,57rO") + 5oXO<5(2A[],7r°'^-AO,7r^'^-4vr°'=) 

Jy '' -V-^ 

=0 by (2.44) 

= [ 5oXO(57rO”. (2.46) 

Jar 

A similar calculation using (2.45) gives the following contribution of terms in¬ 
volving OqX^ in (2.28) after inserting (2.39) there: 


(2) :=-[ dk (doX^dTT^^) =- [ doX^dTT^^ , 

Jy ^ ' Jay 


(2.47) 
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which had to be established. 


□ 


For further reference, we note that those terms in (2.28) which involve diX^ 
and its space-derivatives take the form 


(3) 







+ [ diX^67T^^ . 

Jar 


(2.48) 


We continue with the terms which do not involve X^ and OqX^. Recall 
that Y is given by the equation = 0}. To avoid ambiguities, we will use 
the notation 


Y ;= X^dk , thus X = X^do + T, with dy^{Y) = 0. (2.49) 


It has been shown in [33], for variations such that the image of V in ^ 
remains fixed (equivalently, = 0), that we have 


+dk 



(2.50) 


recall that we allow det gmn to have either sign. Here we use the usual definitions 


pkl 


^ pO _ ^ /|0 

vldit^ , 


(2.51) 

(2.52) 


where g^^ is the n-dimensional inverse of the metric gki induced by g^i, on Y. 
(More precisely, the calculation in [33] has been done in dimension 3-1-1, but 
the same calculation in higher dimensions leads to the formula above.) To avoid 
the overburdening of notation, we will not make a distinction between the fields 
on Y and their pull-backs to V. 

We use (2.50) as follows: Consider any differentiable family of fields t i—)■ 
A{t)). We first note the identity 


dtA%S7r^’'-dt7r^’'6Al, = + 5{7r^’'dtA%) 

+7r^^^[dt,6]Al,. 

Equation (2.50) implies 


dt detgn^n|i^') 

6i7r>^^dtA%) = 

s [^^^^dt{^/\deigmn\K) 


QkidP^'' -I- dk 


IOtt 



gkidtP^^ + dk 


IOtt 



/^rOfe 


(2.53) 


(2.54) 


(2.55) 
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Inserting this into the right-hand side of (2.53), we find 


- dtT^^''5A% = ^ {dtguSP'^^ - d, 


- dtP^^ 


(n-3) 


[dt, 6 ]{y'\det gmn\K) + T^gki[dt, S]P^^ [ 5 ,dk] -rr^^dt 


-[dt,dk] 7 r °°(5 -fin +4 7 r°°[ 5 t,( 5 ] -fin 


-dk dtTT^^S -TiiT - ^ Utt' 


(2.56) 


When dt = dyO = do and all commutators vanish, one obtains 
doAljn^'' - do7r^''6A% = ^ [dogkiSP’^^ - doP^^Sg^i) 


-dk doTT^^d -fin - ^0 -fin . (2.57) 


We note that for any field /, with 5f := dxf\x=o, it holds that 

dtSf = X<^dMx=o = dx{X'^d^f)\x=o-dxX'^\x=o{daf) 

= 6{dtf)-5X^d^f. (2.58) 

While this can be used to analyze the commutator terms in (2.56), it is cal- 
culationally advantageous to proceed as follows: Using (2.39) and (B.6), Ap¬ 
pendix B, and taking an extension of X^ off 'P such that doX'^ = 0 we have 

{^yA%)57t^’^ - = 

= y* - di7:>^''5A%) -dkY>^7T^'^6A% + 2d^Y’^5 J 

'-V-' '---' 

* =2aiF'=5(7r-A0j 

-didkY’^STT^^. (2.59) 

For the term *, we can use (2.56) with dt = dyi and with vanishing commutators; 


d^AljTT^^^-d,7r^^''6A% = ^ (^dmiSP'^'- ^^P'^ 


-dk 9j7r°°(5 


(57r°° . (2.60) 


The middle term is given by (2.50). The term involving can 


be rewritten as 


Jl 4O I iO 4O _ -^Vdet gmn 7 iZ-pO I iO jO 


A'J A'J I A'J '' "'"“a'""- 

^ = TT Aki Ytt Ako = - — - 

\/det gmn ~ii ^ ^ a _00 

=- 16^^’ 


g i ki -hg 


(2.61) 


where we have used the second line of (B.ll), Appendix B: 


— -1 

A-ini^ — “ 


1 4vr°° ^0 9° 


2 vrOO ^‘gOO ' 
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(2.62) 


In the first line of (2.61), and elsewhere, we use the notation 

jyk _g0k/g00 ^ ^ 

Now, we add all terms depending on 
1 

Ihvr 


{y' {dmidP^^ - d. 


- d,.P^^5gki) - dkY^ [6P - gki5P^^) - 2d^Y^5P^ 

= {CygM^fP^^ - CyP^^Sfgki) . (2.63) 

The remaining terms gather to a full divergence. 




-Y^dk 


M 


TT 


00 


9,:7r^^(5 ] — di\ 5Y 


-kO 


* ’ ttOO 


- diY% 


^00 C 

TT 0 


TT 


0k\ 1 


TT 


00 


^iO 


= dk 


) + didkY'^Sin^^N^) 


5i(yv°°) SN’^ - {Y^diN^ - diY^N^) 


(2.64) 


CyN^^ 


leading to 


L , (2,65) 


N 


A case of interest in its own is that of vector helds X which are tangent to 
y. An alternative, standard and rather more straightforward treatment of this 
case is presented in Section 4.2 below. Here we show how the relevant identity 
follows from the calculations above; 

Proposition 2.5 Suppose that the map Tp leaves P invariant: 

'ip{t,P)cP, (2.66) 

and let Y denote the generator of the flow 1 1 -)> Tp{t, •) on P. Then 

{CyguidP'^^ - CyP^^ 


! y 


fdr 


2Y^5P^k - Y'^P'^^ 


- 3271 / X^5£^i 

J-r 


(2.67) 


Proof: We return to (2.27) with A = 0, which we repeat here for the conve¬ 
nience of the reader 


yxAl,)57T^^ - yxTT'^nSA^uu = -2X‘^6£\ + X^£^^^5g 




+Y^9^|<5[v 1^(V^X^ - V^x^)]| + a, \x\>^''5Al, - 




V^(V^(5X^ - V^(5X^) 


( 2 . 68 ) 
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Integrating over 'Y, the last term above leads to the following boundary inte¬ 
grand: 

^ [v¥(V-X0 _ ^ ^ - g-'^g^^^)d,{g,xX^) 

= {2NX^K^i - doX^ + NWiX^ - XWiN^) . (2.69) 

Here, and elsewhere, D denotes the covariant derivative operator of the space- 
metric gij. Now, we have seen that 

[ - {CxTrnSA^J d^A = (1) + (2) + (3) + (4), (2.70) 

where (1) has been defined in (2.46), (2) in (2.47), (3) in (2.48) and (4) in 
(2.65). Let us take an extension X of the vector field Y (which is tangent to Y 
and defined so far only on Y) such that X^ = 0. Then (1) and (3) vanish by 
(2.46) and (2.48). From (2.47) and (2.65) one finds 

[ iCxA%)5f7r^’^ - (CxnnSfA^ 

J-y 

= ( 2 ) + J (y^xQki^fP^^ - YxP^^Sfgki^ 

-J-[ \cx . (2.71) 

lovr Jqx L J 

Using (2.50) and (2.69) to rewrite the right-hand side of (2.68), with some work 
(2.71) can be rewritten as 



The first term in the first line cancels out the first term in the last line, providing 
the required result. □ 


Our considerations so far have taken care of all the terms involving X* and 
doX^. We continue with the analysis of the remaining terms, with the aim of 
proving (2.29). For all practical purposes, the calculations are the same as if 
we wanted to analyze the variational identity (2.27) for vector fields such as 
X* = 0 = doXf^, and so we proceed accordingly. Equation (2.27) specialised to 
this case reads 


{^xA%)67r>^- - {^xxnSA% = 

X^ [(aoA°J<57r^" - (9o7r^")-57l°,] + di{d^X^5x^^^) 

+dkX^ - vrO^F^o)] • (2-73) 


18 



In order to analyze the first term on the right-hand side above we use (2.57). 
Integrating over 'Y gives 

(dogkiSP^^ - 


+ / 


f y 


- 

Jar 


9o7r“5 


TT 


kO 


5ovr°°(5 ( ^ 1 - 5o ( ^ ) dvr' 


TT 


00 


- do 


TT‘ 


kO 


TT 


00 


6tt' 


00 


TT 


00 


TT 


nO 


TT 


00 


,00 


(2.74) 


To analyze the boundary term in (2.74), we start with the identity 


(^) + (^) = (5 


TT 


nO 




yj-OOyj-nnl 


(2.75) 


Writing 


TT 


nO 


r,nO 


q := 


V\ 


TT^^TT^^I \/ \ 9 ^^[ 


and assuming that sgnTr^^ sgnyr^^ = — 1, one finds 

2y^|7r00y^nn| _ — \/\^\\/ g^'^\ = 9a^ 

We will write 


A := kdet gAB 


(2.76) 


(2.77) 


(2.78) 


for the pull-back to dV of the (n — l)-dimensional volume density on the bound¬ 
ary dY. In this notation we have 


^00;: 
TT 0 


TT 


nO 


TT 


00 


+ tt^^6 


TT 


TT" 


nO 


^ (2.79) 


Svr yA+Jf 8 tt 


where a := arcsinh(g) is the hyperbolic angle between the vector orthogonal 
to the hypersurface Tp{{t} x V) and the world-tube {y” = const} (e.g., a = 0 
corresponds to the situation where the vector ip*dt is tangent to the tube). 
When sgn7r°° sgnvr'^” = 1, instead of (2.77) and (2.79) we get 


and 


2v1' 

yj-00yj-nn| 

2 

IGtt 

O 

o 

/vr’^OA 

/ 





TT 


TT" 


nO 


A Sq 


= —Sa , 


Stt y/i — q2 8 tt 


(2.80) 


(2.81) 


where now q = sin a. 

This, and a calculation identical to the one leading to (2.53) imply that the 
boundary integral in (2.74) takes the form: 



doTT^^d 



do 


/^no 
I ^nn 




+ — (9ocu^.^ — 9oA(5a) 
Stt 
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Finally, we obtain 




[dogkiSP'^^ - doP^^5gki) + ^ J^^X\doa6X - doX5a) 


+ / 9fcX° 9o7r“5 ^ - ^0 ^ SP 


+ / XO doTT^^d 

Jd-r 


(2.82) 


We continue with an analysis of the boundary term with X* set to zero: 


X07r^^<5yi;)^ - XV^'SAl^ = X^^^^^''5Al^ . 


(2.83) 


Let us exchange the role of y"' and y^. Identities (2.44) and (2.45) become 
constraints on the boundary of the world-tube T x dV, where T is the time-axis: 


An _ 
■^nn 


— _ : 

na 27r' 


L + <>■*) , 


(2.84) 

(2.85) 


and a, 6 = 0,1, 2,... , n — 1. They imply 


n^^SA:, 

(n - 3) ^ 




5iV\detg,a\L) - —gabSQ^'^ + da tt^^S — , (2.86) 


where Q°'^ and L have been defined in (2.6). Equation (2.86) gives 


(n-3) f 
Ibvr Jq- 
(n-3) f 

167 r .L- 


X^5{y'\detgcd\L)--i^ I X+ / XU7r"",5 ^ 


^0 o nn 


X0<5(| ^/d^\L) [ X^gabSQ^^ + [ X^dA ^,5 ( 

. IdTTjgy Jgy [ \7r^^ J} 


+X/"’ (^^J - So j [.””8o 

We turn attention now to the Komar-type boundary terms in (2.27): 
-^6 [v1^(V"x0 _ _ 1 ^(V"5X0 - V°5X") 


.(2.87) 


= X°5 (7r”^r°o - 7r°^r;)o) + df^X^Sir^^^ . 


It holds that 


TT-^rOo - vrO^F^o = ^V\d^\r^Lao - TT^^do 


( 2 . 88 ) 
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The last boundary term in (2.87) cancels the corresponding term in the 
above formula when we add them together: 


Jar 

(n-3) 


(vr-^rOo - + 


'-v-' 

=(!)+ d'y—terra from (3) , cf. (2.46) and (2.48) 


1 


Ibvr 


tar 




Ibvr 


<ay 


Stt. 


aO 


dr 


+ / + / x^dA 


Jar 

+ / X° 

Jar 


ar 

^nO 


TT 0 


TT 


nA 


TT" 


TT" 


doTT^'S _ - So — Utt^ 


TT 


nO 


7r'‘ 


(2.89) 


Now, a rather lengthy calculation shows that 


1 

—( 
87r 


1 


\/|det^(5“°7vao+-(« - 3)L) 


Ibvr 


dab^Q 


ab 


1 

IGtt 


{21^60, - 2iy^SQA + i^Q^^Sqab) , 


AB 1 


(2.90) 


1 

IGvr. 


ar 


which enables one to rewrite the second line of (2.89) as 

X^{2u5Q, - 2i/^(5Qa + i^Q^^Sqab) ■ 

/ 

Formulae (2.82) and (2.89) give the following intermediate result: 

f X° (dogkiSP'^^ 




+ 2/ X^{doa6X-doXSa) 

Jar 


- [ X^{2v5Q - 2v^5CtA + i^Q^^SgAB) + [ {**) x 
Jar Jr 


= - BaX 

Jar 


TT 0 


TT 


.nA 


TT'' 


+ IGtt / X° {-2d£% + S^^dg^,) , (2.91) 


fy 


where the expression (**) is given by 

M = ^SAl +s{ A>‘t% - - a„ i,««. 

=:(a)/(167r) ■.=^ \det gij\{b) / (IQ-k) '' =:(c)/(167r) " 

(2.92) 

This is precisely the multiplicative factor of d^X^ in the sum of volume terms 
in (2.48) and (2.82). 

Equation (2.91) ends the proof of (2.29) when X^ = const, because then 
dkX^ vanishes, while the first term in the last line of (2.91) integrates to zero. 

To calculate (**), we will need the usual ADM formulae (B.l) for an (n + 1) 
splitting, together with 

r^o = dk log N-—Kik, rgo = do log N + rOk log N - -j^Kik 
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(recall that Kki = -NT^ = ^{DiNk + DkNi - dogu))- Furthermore, 

> ^Oj = D,N^ - NKKj + — (N^Ki, - D,n) . 

Using the above, we find 

(6) = iV/'^r°o - = 2F)^iV + ^ - DiN'^N^ - 

= 2D^N + - 2NiK'^^ , (2.93) 

(c) = go( Vde^gmn )^^fc _ jvfc^(Vd^gn»i) ^ (2.94) 

To calculate (a), the following intermediate results are useful; 

/ \ 

Ago = rgo = - NiD^N\ +ND^N + — i^N^N^Kij - N - N’^DiNj 

N*=+NiD^N*=-2NNiK<'>^ 

(2.95) 

Ak -pfc ^ ^ X^P^ ^ X^pO 

^Oi - ^Oj- 2 ^ OA - -L Oi - 2 ^ o« - 2 ^ oo 

1 / \ 

ylg^. = DjN'^ - - DjN) - -S^ [do log N + DiN^ - NK) 

(2.96) 

/i fc _ -p/c c/c -pA _ t^Ai e/c T^l e/c P'0 

Aj - i ii - (’(ji i)A - i)l - i)0 

= fg, - log iV+—iF,, (2.97) 

We pass now to 

dkX^^ = = dkX^ (7r°°Mgo + 27r°^Mg,. + vr'W^g,.) . 

Using the above, one can hrst notice that (**) does not contain SqN or doN^. 
After some further work we obtain: 

(**) = 5 (f- ,5f,f+ 2D’^NSa/^ + + 

+^DiN{N'^5N^ - N^SN’^) - ^KSN^ 

+ ^ N^6DiN^ - DiN^SN^ - N^S{DiN^) + DiN^dN’^ 

+^nHk - 2N^i^Kf) + . (2.98) 

Those terms in (2.98) which depend on the shift vector but do not involve 
Kij may be simplified as follows; 

-Va^S (fg^. - + 2^d,n{nHn^ - N^6N^)+ 
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-di [Vdetg^n + D^X^dg^^^] +Vdetgmn [g’^^DW.X^ - DiD^X^) dg^i, 

all the terms in {**)DkX^ which involve DN and DDX^, and which we denote 
by (* * *)^ take the form 

(* * *) = DkX^ [N^g^^5 (ff^. - + 2D^N5^ + ^D,N5~g^^' 

= -di [n^ + DkX^dg'^^) 

+2DkX\2D^N5^ + y^DiNSg^^) 

+N^ (^g^^DWiX^ - D^D^X^) Sgu . ( 2 . 102 ) 
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Comparing with the last two lines of (2.101), we see that (* * *) differs from 
the desired expression by a divergence term. Further, all the volume terms in 
(2.29) have at this stage been accounted for. Now, this divergence term leads 
to one more boundary integral, with integrand equal to 

-N^ 

= -N^DuX^ - -g^’^-gi^) 5g,j 

= -2N(^y^DnX^ + DcX^-^yX 

+NX^^^dcX^ (^-^SgAB + V¥^5gn^ . (2.103) 

-VW^9AB&{pr) 

The before-last line above will be part of Cxol. 

Adding all the results above, and using the formula for CxX (see (B.5), 
Appendix B)) we end up with the following identity when X = X^Oq: 


ly 


[CxgkiSP^' - PxP^^ 


+ 2 (yCxoi5X — CxX6a + Xu^dAX^Sa) 
Jar 


dAX^ 


lay 




ygnnj 


N N J 


+ [ X^{2ia5Q - 2ia^6Qa + i^Q^^SgAs) + lOvr [ X^Oa 
Jay Jay 

+167r [ AO {-25£% + £^'^5g^,) . 

Jy 

Using the identity 


TT 0 


TT 


.nA 


TT" 

(2.104) 


N^ = v^ + N^- 


one finds 


dAX^ 


lay 


NX^/^5 


-nA 


+ V^[ 


fN^SN^ N^5N^ 


N 




N 


TT 0 


— 2Az/^(5a 


TT 


.nA 


TT'' 


This, together with some obvious cancelations of boundary terms in (2.104) 
ends the proof of (2.29). □ 


3 The Kerr-de Sitter metrics 

We wish to analyse (2.3) for metrics defined on a half-cylinder [0, oo) x 5^ which 
asymptote to a periodic spacelike hypersurface in a maximal extension of the 
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Kerr-de Sitter metric. (The reader is referred to [9], or [18] and references 
therein, for a discussion of maximal extensions.) Locally, in Boyer-Lindquist 
coordinates [9, p. 102], the metric takes the form 

g = (®dt - + a^)dip)‘^ 

—(dt — asin^(d) d(/ 9 )^ , (3.1) 

where 


P^ 

= -|-cos^( 6 <), 

(3.2) 

Ar 

= (r^ -|- a^) ^1 — “ 2mr , 

(3.3) 

Ae 

= 1-b cos2(6'). 

(3.4) 


a^A 



“ 3 ’ 

(3.5) 


with t G M, r G M, and 6, 9 ? being the standard coordinates parameterizing the 
sphere. 

Above we have used the standard notation for the Boyer-Lindquist 

coordinates. This leads to a conflict of notation, as in all our considerations so 
far t was a parameter along the Hamiltonian flow. The Boyer-Lindquist coor¬ 
dinate t would correspond to the coordinate in our previous considerations 
in cases in which the initial data surface is given by {t = 0}. Note, however, 
that such coordinates can work at most in regions where dt is spacelike, so the 
identification of the Boyer-Lindquist coordinate t with is also to be avoided. 
Our explicit calculations will be done for a boundary surface taken to be 
one of the coordinate-surfaces {t = const r = const}, before passing with the 
boundary to infinity. 

Throughout we will assume A > 0 and a / 0; the case a = 0 has been 
covered in [12], while the case A < 0 is briefly discussed in Appendix E. 

We will keep away from zeros of p, where the geometry is singular, and ignore 
the trivial coordinate singularities sin(0) = 0. Recall that the metric (3.1) can 
be smoothly extended across the zeros of A^, which become Killing horizons in 
the extended spacetime. Note that under our assumptions has at least two 
and up to four distinct real zeros. A projection diagram for a natural maximal 
analytic extension in the case of four distinct real zeros is shown in Figure 3.1, 
see [18] for the remaining cases. The figure illustrates clearly that KdS space- 
times with the corresponding range of parameters contain complete periodic 
spacelike hypersurfaces. Inspection of the remaining diagrams in [18] shows that 
complete periodic spacelike hypersurfaces meeting an infinite number of Killing 
horizons exist only when A^ has four distinct real zeros. In the remaining cases 
there exist asymptotically cylindrical complete spacelike hypersurfaces which 
interect only one or two horizons. 
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Figure 3.1: A projection diagram for the Kerr-Newman - de Sitter metric with 
four distinct zeros of A^, from [18]. The pattern continues indefinitely in all 
directions. 
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When m = 0, an explicit coordinate transformation bringing the metric to 
the usual de Sitter form can be found in [9, p. 102], see (E.2) below, compare [2, 
29,40], 

In what follows the following formula will be useful: 

A := ^JdeigAB = ^e{r‘^ + a?Y - a^/\rsn\^{d) ■ (3.6) 


4 Hamiltonian flows for asymptotically Kerr-de Sit¬ 
ter metrics 


Recall (2.35), which we rewrite as 


/ - C-xP^^^Qk?) + 2 / {CxOdQaSX - CxOd^^XSa) 

Jy ^ ' Jay 


+167r / X^S^^’^dgau - 2Xf^5£' 


ly 


= [ X^{2iy6Q-2iy^5QA + i^Q^^SgAB)+ [ (2Y^6P^k-Y^P^ 

Jay Jay ^ 

=: 9b{S,X). (4.1) 


The subscript ”6” stands for “boundary”. 

Given a compact spacelike surface S and a spacetime vector field X defined 
along S, 9b{S,X) is thus the “non-dynamical” boundary term from our vari¬ 
ational identity. 6h can be thought of as a one-form on the space of vacuum 
initial data. 

As such, we consider the collection of vacuum initial data which approach 
Kerr-de Sitter initial data in an asymptotically periodic end, together with 
first derivatives. (A large class of such initial data has been constructed in 
[15, 20, 23].) Our aim is to calculate 6b for such metrics. 

Now, the value of the integral defining 6b{S, X) approaches the value of the 
same integral calculated in an exact Kerr-de Sitter metric when S is moved 
to infinity along the end. So, the problem of calculating 6b for asymptotically 
Kerr-de Sitter metrics is reduced to one of calculating this integral for exact 
Kerr-de Sitter metrics. 

In spite of various attempts, we have not been able to carry out a direct 
calculation of 6b{S,X) for the metric (3.1) on a general surface t = const, 
r = constk However, one can proceed as follows: 

Equation (4.1) shows that, for all metric variations within the family of 
exact Kerr-de Sitter or Kerr-anti de Sitter metrics, and for all Killing vectors 
X of those metrics we have 


eb{S,X) = Obis', X), (4.2) 

whenever there exists a spacelike or timelike hypersurface P so that dP = 

5U5'. 

In particular, if S' is a sphere 

St,p := {t = T, r = p} 
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in a region where the slices {t = const} are timelike or spacelike, we obtain 

= 0. (4.3) 

Since the Kerr-de Sitter metric is t-independent, it obviously holds that 

dMSt,r,X)) = 0. (4.4) 

We conclude that we can choose any convenient value of t and r to calculate 
db{St,r, X) within each of the regions, where {t, r, 6, ip) form a coordinate system. 
Although this is not immediately apparent, it is the case that the integrand de¬ 
pends smoothly upon the metric when approaching the Killing horizons, hence 
surfaces lying on the boundaries of the relevant regions can also be used to 
calculate 6b{St^r, X). So, in fact, to determine 6b{St^r,X) it suffices to calculate 
the integral at a bifurcation surface = 0. This simplifies the calculations 
enough to make them tractable. 

The reader should keep in mind that when seeking a primitive for 6b on bi¬ 
furcation spheres, variations 5m and 5a have to be accompanied with a variation 
5r of r so that 5Aj. = 0. 


4.1 “Energy” 


As such, the condition = 0 implies that zz 
everywhere in our coordinates. Denoting by 
hnd 


= 0 and Q = 0; recall that = 0 
/a.=o limit lim,^o/A,=,> we 


[ -2u^5Qa = [ -2u^5Q^ = [ 5Q^ = -^^5 (^) , 

JAr=o JAr=o r^ + a^ Ja.=o \ J 

drAr ^ f 47r(r^ -|- a?)' 


lAr = 0 


uQ^^5gAB = 


E{r^ + a^) 

To obtain the above, the following formulae are useful: 

(a^ -I- r^) Ab - Ar 


(a2 _|_ ^2^2 _ ^2 sin2(0)Ar 


—Aj.=0 


2 I 2 ’ 


(4.5) 


^ masin^( 6 *) r, n 9 , n 9,9 9,1 

= ^ 2.4 [(^ - ^ +«)] ’ 




1 _ _ V ArAgp _ 

S\/(a^ + Xe - a? sin^(6»)Ar 


(4.6) 

(4.7) 


i'Q'^^5gAB 


QcDd^^g^^^dAB — —QcD^g^^ 

-v\ {LgAB - Lab) 5g^^ 
uX (^L~g'^'^5gAB + LAB5~g^^) 

b'X (^dr{\fArP) X ^5(A^) -F ^^^^^drgAB5g'^^'^ 
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\P 2/9 


2y/ Aj. 


-driy^)6X, 


Ar=0 2 (^ q 2 _|_ j , 2 ^ 

^/A^^r{\/A^) = -drAr = r -m - 5 -( 2 r^ + a^), 


(4.8) 

(4.9) 


_ sin(0) 
We also note 


2 3 

A0{r‘^ + a2)2 - a2ArSin^(6') =Ar=o (r^ + a^). (4.10) 


1 


1 


ma 


o_ / ^ V ~ O- I Q'y’ ~ ^^2 ’ 


Stt 


/S2 


Stt 


'S2 


(4.11) 


where the last equality can be used to simplify the calculation of the term 

[ <5Q^. 

Ja. 


lAr = 0 


From 6Ar = 0 we get hdrArSr = rSm — (1 — ^r^)a5a. Finally, we are led 


to 


L 


Ar = 0 


2i'5Q — 2 u^SQa + I'Q^^SgAB 

drAr 


1C ® i^r^r r 

= lOTT—-^ 0 ^^ + 47r - ^0 


+ a 

IOtt / m 

lO 


2 '=r2 


"(r^ + a^) 




^3/2 


(4.12) 


The last equality is not completely obvious, and its proof proceeds as follows: 


a ma 

r^T^ + T 


1 drAr /^2 _L „2 


+ a 


2 ^^2 


4 H(r2 + a^) 


a ma 1 
ro-=w + tt: 


+ a 


2 ^^2 


2 + a^) 


+ a 


{r6r + a5a) + 


1 rA 1 

= --g^ada + ^ 


m 


= 5( ^ ) + I r - -r - 


A 


2ma^ 1 
' 2 

2mr^ 1 \ 1 „ 1 


4 H 

^^.(2r^ ' ^2, 
6 


2 + -(r - m) - —(2r^ + a^) 


A 


r2 + a2 2 




(4.13) 


r Ar 

r^+a^ 

where the following identity is useful: 

—^a<5a = -S- . 

3H2 2 S 

Equation (4.12) is equivalent to the statement: 
Proposition 4.1 The flow generated by the vector field 

X = VEdt 

is Hamiltonian, with Hamiltonian Hx given by 

® »vsa, = A 75 ■ 


(4.14) 
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4.2 Angular momentum 

There is a standard calculation which shows that 6b{St,r,d,p) is, up to a multi¬ 
plicative factor, the variation of total angular momentum, which we reproduce 
here for completeness; compare Proposition 2.5. For this, let X be any vector 
field tangent to V. The divergence theorem and the vacuum vector constraint 
equation = 0 give^ 

-J:=- [ P^jX^dS, = - [ P^WiXj = ^[ P^^^x9ij . (4.15) 

iJa-y iJ-f 

Hence, if X does not depend upon the metric and is further tangent to dP, 

-5J = ^[ 5P^^^x9ij + P^^^x69ij 

Jy 

= ^ [ SP^^^x9ij+-^x{P^^6gij) -^xP^^Sgi, 

Jy '' -^ 

Dkixkpasgij) 

= ^[ 5P^^^x9ij--^xP^^dgij) 

^iJy 

+7^ [ P'^Sgij X^dSk , (4.16) 

^7 Jar '—' 

0 

providing the usual Hamiltonian formula for angular momentum: 


-6J 


1 



Spij 


dpij 

dip 


^9ij • 


Comparing with (4.1), we obtain an alternative proof of Proposition 2.5 for 
vector fields tangent to S: 


db{St,r,d^) = -2^6J . 


(4.17) 


When X remains tangent to P but is allowed to depend on the metric, there 
arises a supplementary term 

^ [ P^^^sx9ij = - [ P^jdXUSi. (4.18) 

277r iJay 

We see that whenever 5X is a Killing vector (which is the case in our consider¬ 
ations here), or when 5X vanishes at dP, the integral (4.18) vanishes. 

From (4.11) we obtain 


/ 771 (1 \ 

9b{St,r,d^) = -2j5J = -16x6 ■ (4.19) 

^The minus sign in (4.15) has been intriguing to us. In this context it is helpful to realize 
that in a field theory where the energy density is given by Too, energy-conservation in its usual 
form requires the momentum density to be defined as —Toi- A simple consistency check is 
provided by a massless scalar field of the form 0 = f{x — vt), where the minus sign is needed 
to have the field momentum positively directed along the velocity vector. Compare [7]. 
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4.3 The Henneaux-Teitelboim vector field 


In [29], in the closely related case of negative cosmological constant, the authors 
consider the vector field 

Xht '■= dt + ■ (4.20) 

The volume part of the variational formula (2.3) is linear in X, which gives 

0b{St,r,XHT) = 0b{St,r, dt) + -^0b{St,r , dip) . (4-21) 

Hence, 


db{St,r,XHT) 


IGtt 





(4.22) 


We see that Xht is not Hamiltonian. 

On the other hand, (4.22) and our remaining considerations so far show 
that: 


Proposition 4.2 The vector field 

X = '^Xht = Sdt + -^dp (4.23) 

generates a Hamiltonian flow on the space of metrics, with Hamiltonian Hx 
given by 

TYl 

= (4.24) 


4.4 Kerr-Schild coordinates 


In [9, p.l02] it is emphasized that the Kerr-de Sitter metric is asymptotically de 
Sitter as it tends to the de Sitter metric in the limit as r goes to infinity. One way 
of making this precise is to invoke the Kerr-Schild coordinates, with a de Sitter 
background metric. Following [2, 26] we use the coordinate transformation 


dr = 


~dt -\- 


2mr 


l-r^)A, 


dr. 


2mr a 


to obtain 


with 


aA 

= ''■f - + (,.2 + „ 2 ) A , 


2mr,, , ...o 

9 = 9ds^ - —{kpdx>^) , 

P 


dr 


(4.25) 

(4.26) 


9dS = 


1-2^ A 


-dr^ + 


(l - (r2 + a2) 


dr^ + 

Xe 
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(4.27) 


I (r^ + Q^)sin^(g) ^^2 

k dx^^ = ^dr + -, - ^ - dr - . (4.28) 

^ (l - r^) (r2 + a2) 

The metric gds is the de Sitter metric in unusual coordinates, which can be 
verified by using [9, 26] to find the transformation 

2 r^Ae + sin^(6') 

i?^sin^(0) = sin^(6>), 

T = T, 

4> = 0 


between (4.28) and the de Sitter metric in static coordinates: 

A 1 

9dS = -(1 - -^)dT^ +- xi^dR^ + (d02 + sin2(0)d<h2) . (4.29) 

3 1- 

The vector k^dfj, is null both for g and gds- 

According to the above, the natural “time-evolution” Killing vector Xks in 
the Kerr-Schild setting is 

Xks ^dr = Edt + . (4.30) 

This coincides with the vector field (4.23), with the associated Hamiltonian 
given by (4.24). We can also calculate directly: 

db{St,r, Xks) = 




Eeb{St,r,dt) + —9b{St,r,d^) 


IGtt 




m 

■=3/2 


Aa . / ma 

- t" 


IGTrSi^ 


(4.31) 


5 All Killing vector fields generating Hamiltonian 
fiows 

We have seen that the “energy functional” 



generates the flow of the vector field V^dt, while the “angular momentum 
functional” 



generates the flow of —d^p. This allows us to prove the following: 
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Theorem 5.1 A vector field X defined along an initial data hypersurface dX 
whieh asymptotes to a Kerr-de Sitter Killing veetor field as one reeedes to in¬ 
finity along an asymptotically KdS end, 

X ^ -\-fid^p , (5-1) 


where x o.nd /3 are allowed to depend upon J and E, hut not upon coordinates, 
generates a Hamiltonian flow on the space of initial data sets with one asymp¬ 
totically Kerr-de Sitter end if and only if 


^ ^ = n 

dJ dE 


(5.2) 


If (5.2) holds, the Hamiltonian Hx = Hx{E,J) generating the flow of initial 
data associated to X is the unique, up to a constant, solution of the equations 


OeHx = X, djHx = -fi . 


(5.3) 


Proof: Let us denote by Tr C ^ a family of domains with smooth bound¬ 
aries Sr such that Sr is homologous to spheres of constant t and r in the 
asymptotically KdS region of , with UrTr = y. 

Let Xi denote any vector field which asymptotes to , and let X 2 denote 
any vector field which asymptotes to dp. Then the vector field 

X 3 := X - xXi - /3X2 


asymptotes to zero. 

We have seen that, for all vacuum variations of asymptotically KdS initial 
data it holds that 


lim X [ (CxigkidP’^^ - £xiP^''5gki) + 2 lim x / {PxOdgOidX - Afe) 
R^OO V J R^oo Jsr ^ ^ 


= lim xOb{SR,Xi) = xdE , 

R—>-oo 


(5.4) 


lim/?/ (Cx 29 kiSP^^ - Cx 2 P’'^Sgki)+2 lim f3 f {CxOgadX - CxOgXda) 

R-^00 V J R-^00 2 u 2 u 


= lim fl9b{SR,X2) = —fldJ. 
R^oo 


(5.5) 


lim 

R^oo 




Cx.gkidP^' - Cx,P^' 




'Sr 


= lim 9b{SR,Xs) = 0. 
R^oo 


(5.6) 


(In the first two equations, the boundary integrals involving a and A vanish 
because dt and dp are asymptotic Killing vectors; in the last equation, all 
boundary integrals vanish because X 3 asymptotes to zero.) It easily follows 
that 

xSE-fl6J = (^CxgkiSP^'- PxP^^Sgki) (5.7) 
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The left-hand side will be equal to SHx, for some function Hx = HxiE, J), if 
and only if (5.2) holds. □ 

The theorem implies that neither the field dt, nor the Henneaux-Teitelboim 
field XhTi generate Hamiltonian flows on the space of metrics. 

A corollary of Theorem 5.1 is that every Killing vector field of the Kerr-de 
Sitter metric can be rescaled by an m- and a-dependent factor so that the 
rescaled field is Hamiltonian. Indeed, consider a vector field of the form (5.1) 
which does not satisfy (5.2). Then the vector field 

^{xV^d/dt X j5d/d^p) (5.8) 


will be Hamiltonian if and only if 7 is a solution of the ODE 


d ^ d dx 9/3 \ 


(5.9) 


This equation can always be solved locally using the method of characteristics. 

Note that many such rescalings exist; Suppose that X has been rescaled 
so that 7 X has a Hamiltonian H.yx- Then for any function / the vector field 
f {H^x)lX is Hamiltonian, with Hamiltonian 


Hf'^x = f{Hyx) ■ 


As an example, we note that the vector field 


dt -^9(^=“-Tks with a replaced by its negative” 

3.^ 


is not Hamiltonian. But the Killing vector field 


V2E - 1 

with Hamiltonian —\/2S — 1. Similarly ( 9t — '^9,^ ) is, with 


dt - —9^ I IS, 


Hamiltonian 


m 


(2H-1). 

Further examples of Hamiltonian vector fields include 

da da J 1 

-^Xks - - -J^Xks - . 


generated by H = a = J/M, and 
dm dm „ 


1-^Uks-^9, 


M2 


3M 


(5.10) 


(5.11) 


generated by H = m = M[1 + AJ 2 /( 3 M 2 )] 2 . 


6 Black hole thermodynamics 

“Black hole thermodynamics” can be thought of as considerations based on 
variational identities involving global invariants such as mass and angular mo¬ 
mentum and local invariants associated to event horizons. In stationary space- 
times with a preferred “stationary” Killing vector, say X, and a Killing horizon. 
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one proceeds as follows; One defines the total energy E as the Hamiltonian as¬ 
sociated to the flow generated by X. If X is tangent to the generators of a 
horizon with compact cross-sections, a rather general calculation shows that 

6E = ^SA, (6.1) 

OTT 

where A is the area of the section of the horizon, and k is the surface gravity. 
If, however, X is not tangent to the generators, then there must exist another 
Killing vector X tangent to the generators. In spacetime dimension four, mild 
supplementary conditions guarantee that one can find a linear combination, say 
r] = d^p, oi X and X which has 2it periodic orbits: 

X = X-nr]. (6.2) 

Note that a specihc normalisation of X, determined by X, has been chosen in 
(6.2). Equation (6.1) is now replaced by 

dE = ^5A + n5J . (6.3) 

OTT 

As already discussed in e.g. [25], a basic problem with this program for 
Kerr-de Sitter spacetimes is the choice of the Killing vector X above. Let us 
illustrate this with some examples. 

6.1 dt 

Suppose, first, that we decide to use the 9t-vector held of the Boyer-Lindquist 
coordinates as the preferred Killing vector held X. As we have seen that this 
Killing vector is not Hamiltonian, the prescription cannot be applied. This last 
problem can be hxed by taking instead X = V^dt, which is associated with the 
Hamiltonian We will return to this choice in Section 6.3. 


6.2 dr 


Yet another choice, which is directly Hamiltonian without the need to do any 
rescalings, is the Kerr-Schild vector held Aks of (4.23), and let us therefore 
investigate this choice in detail. 

First, some general considerations are in order. In maximally extended 
Kerr-de Sitter spacetimes, the Killing horizons are located at 


= (r^ -I- a^) 



— 2mr = 0 . 


Let r* be a root of this equation, then it is easily seen from (3.1) that the Killing 
vector held 


A* — ^ dt 


O O 


_ ^ , (2mar*) ^ 

- Aks + ^ Op 


dp = Aks + 


A 


+ ri 


— a- 


dp = dr + 


g (3 - Ar^) 
3 (a2 -h rl) 


du: 


= : 0 * 


{o? + rl) 


(6.4) 
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is tangent to the generators of the corresponding Killing horizon. (In the last 
equality above the equation satished by r* has been used.) 

We have seen that the Hamiltonian mass associated with Xks is 

777 

Mh = ^, (6.5) 

and that the total angular momentum J is 

777 n 

j='^ = MHa. ( 6 . 6 ) 

Using (3.6), the area of the horizon is 


= 2tt / A|r=r* d9 = 


27r(r2 + a^) 


(6,7) 


As such, the surface gravity of a Killing horizon can be dehned 
through the formula 

df,{X^Xa)U,=-2K,X^. ( 6 . 8 ) 

A convenient procedure to calculate k* proceeds as follows: let b be any one- 
form which extends smoothly across the horizon and such that b{X) = 1. Then 
K* can be obtained as the value at the horizon of minus one half of the left-hand 
side of (6.8): 

I-'-' — - (6.9) 


= -- 6 (V(X“A„))|^, . 

Note that the left-hand side is independent of the choice of b, and so is therefore 
the right-hand side. 

Recall [9] that an extension of the Kerr-de Sitter metric (3.1) across a Killing 
horizon can be obtained by introducing new coordinates 


dv = dt + [ 


+ 0 ? 
Xr 


dr , 


d(p = dcp + E—dr . 


( 6 . 10 ) 

( 6 . 11 ) 


Hence the one form 


b = —dt + 


Ar 


dr 


extends smoothly across the Killing horizon and satishes 6(A*) = 1. Thus 

(r^ -|- a^) 


= — lim 

= — lim 


2 A^ 
(r^ -|- a^) 


r-^r* 2Ar 
drAr 




g''''drigiX,,X,)) 


Here we have used 


2(r2 -ha2) 

A. 


( 6 . 12 ) 


9 = 


r2 -|- a2 cos2 9 ’ 


A^ (r^-|-cos^ 0) 2 n 

ff(W,W) =--^-^ + 0 ((r-r,) 2 ) 

(a 2 -F r 2 ) 


(6.13) 
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We claim that we have, consistently with what has been said above, 


A . 
1 - r‘ 


4 (j.2 _|_ ^,2^ 

5J 




m 


^2 

SMi 


= H^), (6.14) 


O. k*/ 2 A*/47r 

where r = r^. Indeed, our variational identities lead to (compare [ 8 , 27, 36-38]) 


SMh = — k^6A + . 

87r 

To prove (6.14), we rewrite (4.12) in the form 


aV.=. ^'ma 1 

-o ^ __lO H“-^ 

+ r? 4 VH(r2 + a^) 




= 5 


m 


53/2/ • 


(6.15) 


(6.16) 


Dividing by v S, we are lead to 
a ma 1 

0 ^^ + o - ^0 


4E{r^ + a‘^) 

Inserting the identity 




= 5 


m \ m 

W) + 


m Am A 

J 


9 Cl ^ TTICI 9 

r 


drAr 


4H(r2 + a^) 


+ a? 


(6.17) 


(6.18) 


+ a? 

into (6.17) one obtains (6.14). 

6.3 ^/Zdt 

The dynamics associated with the vector field \/S9t of Proposition 4.1 is gen¬ 
erated by the Hamiltonian 

Mh := mH-3/2 ^ 

We have the decomposition 


V^dt = Vh ( + 


“ .ad - “ 


a^ + rf'* a2 + r2 ” 


=:r 2 » 


5^ = , (6.19) 


= :X* 


where W is tangent to the generators of The surface gravity of A* is 
with K* given by (6.12). By inspection of (6.16) one obtains; 


5Mh = — k^fSA + . 

87r 


( 6 . 20 ) 
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6.4 

We have seen so far two sets of variational identities of thermodynamical type 
which are satisfied in Kerr-de Sitter spacetimes, and one can clearly produce 
an infinite collection of such identities by considering all possible Hamiltonian 
dynamical systems associated with all Hamiltonian asymptotic Killing vectors 
of Theorem 5.1. It is thus natural to raise the question, whether any such 
identities are singled-out by the geometry. It turns out that this is indeed the 
case, and can be seen as follows: 

Recall that the equation = 0 can have up to four distinct roots, and that 
there are always at least two distinct zeros when no naked singularities occur. 
Further, there is always a negative simple root. Consider, then, a Kerr-de 
Sitter spacetime with a a^A < 3, and choose a non-degenerate Killing horizon 

= {r = r*}. There exists precisely one preferred Killing vector associated 
with this horizon, namely the unique Killing vector which is tangent to and 
has surface gravity equal to one: 


X.-j^ .— X.-i^ 


2{rt + a^) 

drAr 





a 

+ 



( 6 . 21 ) 


where A* has been given in (6.4), and we have used (6.12). 

It turns out that A* is Hamiltonian. To see this, we rewrite (6.15) as 
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— {SMh — H*(5J) . 


( 6 . 22 ) 


This means that the Hamiltonian A^/Stt generates the vector field 

— (-Aks + = K~^X^ . 

Equivalently, the flow of is generated by 

A = A±A . (6,23) 

Stt 2“ ^ ^ 

Consider, then, two Killing horizons with areas Aa and generating Killing 
vectors A^. Suppose that Ai is non-degenerate and normalise Ai to unit surface 
gravity. We will denote by H 21 the relative angular velocity of the second 
horizon with respect to the first one: 

Ai = 

A 2 = 

Ai = 


drArlj.^. 


-H'l 


-^r\r=ri 

—-((ri + a‘^)dt + ad^) , 

Or^r\r=r2 


uri-^r\r=r2 

(rf + 2aS(r| - rf) 

{rl+a^)drAr\^^^^ ^ (’'2 + a^)^r’^»’lr=ri 

'- . -' '- .. — 


= :/t21 


—Q 21 


(6.24) 
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We then have the variational identity 


<5 



H21 

Stt 

K.2 


5A2 + ^ 21 ^ J 

9.2-^ 


Sttki 


5 A 2 + 


Ki 


Uj , 


(6.25) 


with Ki given by (6.12), and 11* given by the last line of (6.4). 

For Schwarzschild-de Sitter spacetime, the areas A^ are determined by rj, 
and we note the following relation between the roots of the polynomial Ar|a=o 


r\ + rir2 + = ^=: 

which allows one to express explicitly Ai as a function of A 2 . 

In the general case a / 0, (6.7) can be rewritten as 

I j2 

A^ = * --r ; recall that M = and J = Ma, 

This can be solved for r* as a function of , M and J: 

2 _ {fM'^ + J2)A - 47r£2j2 
~ 47r£2M2 ■ 

Rewriting the equation A,, = 0 as 

(r^ + a^)2 — Arn^r"^ = 0 

one obtains, dropping stars; 


(6.26) 


3A-h 

(6.27) 

(6.28) 


(6.29) 


z^{l - zf - Azu^ - Af{z - 1) = 0, (6.30) 


where 

2 : := A/{A'n^'^), ^ := i~^M, and j = £~‘^J. (6.31) 

Let Za, a = 1,2, be two distinct roots of (6.30). Eliminating /j, between the 
resulting equations, one finds the relation 


Z1Z2 (1 - 221 - 2 z 2 + Z1Z2 + zj + zl) = Aj"^ . (6.32) 


This can be explicitly solved for 22 (-^i), but the resulting formulae are not very 
illuminating. A detailed analysis of the resulting equations will be presented 
elsewhere. 


A Notations 

We summarize some of our notations, which largely follow [33] with, however, 
some exceptions; 

The coordinates are always local coordinates on the spacetime ^. In 
Section 2 the coordinates {t,x) are local coordinates on M x S. However, t 
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denotes a Boyer-Lindquist time coordinate when discussing Kerr-de Sitter and 
asymptotically Kerr-de Sitter metrics. 

The coordinate is constant on 'f' := iI>q{V). 

The symbol , with i,j G {1,..., n}, denotes the n-dimensional inverse of 
the metric gij induced on the level sets of Strictly speaking, we have a 

one-parameter family of metrics gij{t), but we will not use this notation. We 
use the symbol D to denote the covariant derivative operator of gij. 

The coordinate y” is constant on di^. 

The symbol with a,b ^ {0,1,...,n — 1}, denotes the n-dimensional 
inverse of the metric gab induced on the level sets of y^. 

The symbol g^^, with A, B £ {1,... , n —1}, denotes the (n —l)-dimensional 
inverse of the metric gAB induced on dY. 

B Lie derivatives of geometric fields on Y 

The aim of this appendix is to derive the derivative operators which are obtained 
by restricting the Lie derivative of spacetime objects to a non-characteristic 
hypersurface. For composite objects such that a or A, the method is to use 
systematically the chain-rule together with the spacetime expression for the 
Lie derivative of the relevant components of the metric. We use the ADM 
decomposition of the metric; 

ffofc = iVfc, iV2 = yoo = - iV2. (B.l) 

Here, as elsewhere, we let X = + Y, with Y = X^di. We will use the 

symbol Cx to denote the restriction to the hypersurface {y° = const} of the 
spacetime Lie derivative operator, while is the usual Lie derivative operator 
on the level sets of the function y^, viewed as an n-dimensional manifold of its 
own. 

B.l The induced metric 

We start with the following straightforward formulae: 

Cx9ki = X^dogki + NkdiX^ + NidkX^ + Cygki, (B.2) 

= X^doY - Yd^X^ - Yd^X° + CyY- (B.3) 

B.2 The volume density A 

We pass now to 


A^ = dety^s = dety^,, , 


which can equivalently be written as 






N 
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The definition of the Lie derivative 

= X^dxg^’^ - 

together with 

^x\l\ dets-a/sl = (B.4) 

gives the following for X of the form X = X^do: 

CxX = X^doX + = X^doX + Xu^OaX^ . (B.5) 

For further reference we have 

= dx (^X\>^''^ - . (B.6) 


B.3 The angle a 

Using the definition (2.76) of q and (B.3), we get 

Icxq = - —^xg^^ , 

q ^ gOn ^ gOO ^ gUn ^ ’ 

which gives, again for X = X^d^, 

Cxq = X^doq - 


~nk 


vlFV 

and finally for q = sinh a we have the following: 


fdkX^ , 


Cxoi = X^doa — 


g^^dkX^ 


coshay'lg®^ 


(B.7) 


vw^\ 


(B.8) 


B.4 The ADM momentum 

The definition (2.52) of can be rewritten as follows: 

:= x/d^^ - K’^^) 

1 


= X 


(B.9) 


~mk~nl 


~mn~kl 


(summation over spatial indices m, n only), where we have defined the purely 
spatial tensor density 


Xmn ■— I det g^g \ ^mn — \/gmnX-n 


Now, when X has no zeros, P^^ can be defined as a derivative with respect 
to time, as calculated in a coordinate system in which X = dt- This definition 
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implies the Leibniz rule. Applying it to (B.9), we express in terms of “dots” 
acting on and i.e. in terms of Lie derivatives of the metric and of the 
connection. To obtain the general formula, valid for any vector field, we use 
( 2 . 22 ) to obtain 




/il/ 




= + Ttd,x^ + vPdaX^ - k,b„x^ + x^d^v, 


CxTl = dudiX^ + Tl^diX- + Tl^diX^ + TlduX- + TldkX 

-ThdoX^-TlidnX^ + X^dA- 

Recall that = L^q — Using the constraints (2.44) 

+ 2A2,7r°'^ 


(B.IO) 


1 


4° = - 

^ofc 27r00 


1 dk (g°°V|detga^|) _ .po ^ 

2 -00 /I U-+ - „l nOO 


c/OO^I det gap 




1 4 vTdet 5 ^ _ 1 gfcffQP _ ^0 g' 

2 v^ldet gag\ 2 


01 


00 


kl gOO ’ 


(B.ll) 


together with the formula 


pA _ 

fcA — 


dk\/\ detgqgl 
vTdetft^ 


(B.12) 


we find 

The identity 
implies 


pO _ 1 9kg' 

k ko — ~~ 


00 


2 g' 


00 


-r 


„0l 

0 ^ 
kl gOO • 


gmO^kl gmrP'kl — ^nkl — gmnXkl 


nO 

■pn _ T^n r, pO _ ■p^ _L _ _tO 

^ kl ^ kl 9 9m0^ kl ^ kl ' ^00 ^ ’ 


(B.13) 


where T^^ describes the spatial connection on hypersurfaces {t = const}. In¬ 
serting (B.13) and (B.12) into (B.IO) we obtain 




= DmDnX^ + x^dA^ + rlAx'^ + ro - ro 


-?00 KndkX^ + rlkdnX^ + rlkdmX^) 


Ok 


{dmg^^dnX^ + dng^^dmX^) , 


2g' 


00 


where D denotes the covariant derivative in each hypersurface {y^ = const.} 
separately. Using now (B.4), we are led to 


CxIC 


mn 
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= X^d^lCmn + \J\ detgagl DmDnX^ 

+dk (^X^X mn^ H“ J^mkdn^ ~\~ 

O/c 

{^mndkX^ + XmkdnX^ + Xnkd„,X^) 

(a™<,“a„x» + a„s»»a„.jf»). 

Using ADM notation, we have ^y\ detgag\ = A" ^/\detg^\ and 

^xK^mn = X^doXmn + Ay^j det S'mnl DmDnX^ + Cj^Xmn 
+N^ {XmndkX^ + XmkdnX^ + XnkdmX^) 

+ \/| detgmn\ {dmXdnX^ + dnNdmX^) . 

On the other hand, 

r ( 7,^1 7mn ~kl \ ~nl p ~mk , ~mk p ~nl ~kl p ~mn ~mn p ~kl 

Xx \^g g - g g j = g Lxg + g Lxg - g Lxg - g Lxg ■ 
Consequently 

cx {g^^r' - = x^do [r’^r' - r^-g^') + - r^g^') 

_~mn (^pjk~lr ^ _ gkl ^pjm~nr ^ pjn~mv^'\^ g^j^O ^ 

Finally, we obtain; 

CxP’^' = Cx{}Cmn{r’'9^'-r^~g’^')} 

= + Ay^drt^^ 

+A^ (^P'^^drX^ + X\D^X^ + X\D^X^ - 25^'/C™9^ao) 
+\/|det5mn| {d^ND^X^ + D^ND^X^ - 2f^D'^Nd^X^'^ 

- {k.\ + A'g’"") + /C'^ [X^f^ + X^g'^^'^ 

-X\ {X^g^^ + A'5^") - (2A”^/C"^)} drX^ . 

After obvious cancelations we obtain the desired formula: 

CxP^^ = X^doP^^ + C^P^^ + AVI det gmn\ [d^D^X^ - 
+ {vldit^ (^^D^X + -g’^^D^X - 2~g’^^D^X^ 

+ pkl^r _ pkrj^l _ plr^k^ g^^O ^ 
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C Stationarity with respect to variations of the field 

X 

One of the consequences of Theorem 2.2 is, that variations of X enter formula 
(2.3) with vanishing coefficients. In this appendix we give an alternative proof 
of this, under the supplementary assumptions that X is a space-time vector 
field, which is everywhere transverse to X, with v without zeros, and assuming 
variations of the map V’ which leave Y invariant. 

Given a field X on ^ which is transversal with respect to a hypersurface 
Y, and give a variation field 5X, we choose an arbitrary 1-parameter family 
X(A) of (transversal) vector fields such that X(0) = X, fulfilling: 

X(A). (C.l) 

A=0 

A possible contribution of the variation of X to the right-hand side of formula 
(2.3) is linear with respect to 5X. We denote it by A{5X). 

Let be the local flow generated by X(A). We set 

M X X ^ x) = . (C.2) 

Next, given a field configuration g on ^, we define a 1-parameter family g(X) 
of metric tensors as 

5(A) := i(p*x)~^g </>a5 (A) = (p*og{0) = g. (C.3) 

This means that in local coordinates (t, x*) in a neighborhood of Y as in (C.2), 
the metric coefficients g^uW of g{X) do not depend upon A: 

g/Mui^) = g/Mu ■ 

So, from the point of view of the manifold 

MxSdMxGpsMxIL 

and the above coordinate system on it, the resulting variations satisfy 

g^iyi/^) — 0) 

A=o 

and, therefore, all the terms in the right-hand side of formula (2.3) vanish. 

Now, it has been shown in [33] that (2.3) is coordinate invariant. This 
implies that the right-hand side of (2.3) vanishes when calculated in any co¬ 
ordinates (x“) on not necessarily the ones adapted to the flow of X(A) as 
above. Observe that by (C.3) we have 

9 = -^zg, 

A=o 


A ^ 


A=n 




^g/iu — ^ipgfii/ — 


d 

dX 
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where the Lie derivative .^zg of the metric g is calculated with respect to the 


field 


x) 


d 


A=0 


Hence, the sum of all contributions of to the right-hand side of (2.3) is 

canceled by the contribution A{5X) of SX. In other words; A{5X) is equal to 
the contribution oi 5g = SAzg- 

In fact, this last contribution vanishes identically. To show this, observe 
that Z vanishes identically on X because 4>\{0,x) = x for every x £ X. This 
implies that variations of all the Cauchy data, i.e. Sgu, and, therefore, also 
5\ and 6gAB, together with Cauchy data of matter fields: 5ip and 5p, vanish 
identically. The only non-vanishing contribution could, therefore, come from 
6Q, 6Qa and da. Hence, the total contribution of these fields, according to 
(2.3), is equal to: 


^((5X) = ^J (^X6a + Z/5Q - • (C.4) 

These quantities obey, however, the following identities: 


<5Q = 

(^X-dAiXB^ij da, 

(C.5) 

+ P\ = 

—XdAa , 

(C.6) 

dQA = 

—XdAda 

(C.7) 


(see Equations (7.11)-(7.13) in [33]). Inserting (C.5) and (C.7) into (C.4) and 
integrating by parts they cancel each other, so that A{6X) = 0. 

For further reference, we note the following formulae for the variation of the 
remaining fields involved; 


dgAB = 0, 

6goo = 5(X|X) = 2{X\5X) = -2v6v + 2ua5u^ , 
^goA = d{X\dA) = gABdv^ ■ 


Moreover, using the invariance of Q = and Qa = <5%, we have: 
0 = ,5 (z/QOO) = + vdQ^^ , 

0 = dQ\ = d {Q^^goA + Q^^gBA) , 


whence, 

lA 

^ _~fA ^ ^ gOO (^^u^du - dB^^ . 
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D Spacelike vectors in adapted coordinates 

In this appendix we show that the variational formula (2.3) for vector fields X 
tangent to the initial data surface coincides with (2.34). For simplicity, and 
consistency with [33] we assume variations that satisfy the linearized constraint 
equations. 

Using (2.90), one can define a Hamiltonian by performing a Legendre trans¬ 
formation which leads to the so-called “purely metric” formula (9.1) in Ki- 
jowski’s paper [33]: 

+ ^ / (D.l) 

lovr Jqy 

where 

^=-^1 Q^ = -^[ fV|det5cd|(r°^aO+^(n-3)L)l . (D.2) 

oTT Jqv ovr Jqy I 2 

For the purpose of the calculation here, it is convenient to change the nota¬ 
tion so far as follows: Instead of we write Py , and instead of <5“^ we write 
P^. This notation emphasises the fact that those fields describe the “ADM 
momentum” of the surface V, respectively the “world tube” T obtained by 
flowing dV along the vector field X: 

r := {Q^{x)\x G dV}. 

In the case of current interest, where X is tangent to V both surfaces coincide, 
hence so do the corresponding ADM momenta. 

In this notation, (D.l) takes the form 

h^L Iv i 

So far the variation SX of X was assumed to be zero, and we used adapted 
coordinates so that X = do- The reader is warned that do has nothing to do 
with a time coordinate in space-time, but is related to a parameter along the 
flow of X. Now, we want to rewrite the formula in a way which allows variations 
of X. For this purpose observe that the upper index 0 in (D.3) describes the 
transversal, or “normal” (with respect to dV) direction, which we denote by n, 
whereas the lower index 0 describes the direction of the field X. Hence: 

Pt\ = Pt\XK 

This leads to an associated rewriting of (D.3): 

xHPt\ = I (P^'Sgki-gkiSP^^)+^ [ {X6a - a6X) 
JdV Jv ^ / OTT Jqy 

+ [ X'^pJfdgu. (D.4) 

Ibvr Jqy 
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Observe, now, that for X tangent to V we have a = 0 (see formula (2.5)). 
Using further = Py, we are led to 

^ 2XHPv\ - X^P^^Sgu = ^ (Pv^gki - guSP^') , (D.5) 

where a “dot” denotes the Lie derivative with respect to X. This is a partic¬ 
ular case, with V = T, of formula (2.3), and coincides with (2.34) within the 
collection of solutions of field equations, as desired. 

For further reference, we note that in fact we also have the pointwise iden¬ 
tity: 

2,ySQ-2iy^SQA + !yQ^^SgAB = 2 XHPv\ - X^P^^5gki. (D.6) 

which is equivalent to (2.90) when X is tangent to P. 

E Negative A 

There exists a clear prescription how to calculate the Hamiltonian mass of a 
family of metrics asymptotic to a fixed background metric [10]. This is the case 
for Kerr-anti de Sitter metrics with negative cosmological constant, which are 
all asymptotic to the anti de Sitter metric [13,14,21,29]. We emphasise that 
this the key difference between A < 0 and A > 0, as considered in this work: 
there is no single metric to which the Kerr de Sitter metrics converge as one 
recedes to infinity along asymptotically periodic ends of initial data sets. 

More precisely, we use the standard form of the background anti de Sitter 
metric, 

b = -(^l + 

where, as usual, = —3/A. We consider the space of initial data sets for the 
vacuum Einstein equations which along P’ := {T = 0} approach the initial data 
for b as R tends to infinity at a rate made precise in (E.ll) below. We wish to 
check whether the Kerr-de Sitter metrics can be put in the relevant form, and 
calculate their Hamiltonian mass. 

In this appendix we apply this prescription to the Kerr-anti de Sitter met¬ 
rics. As such, for such metrics the mass is not a global invariant anymore, but 
the component of a linear functional on the set of KIDs for the anti-de Sit¬ 
ter metrics [11,17], which transforms as a Lorentz covector under asymptotic 
isometries of the anti de Siter background. We will ignore the remaining com¬ 
ponents of the functional and consider only the “energy component”, since the 
transformation properties of the associated object are well understood. 

When A < 0 and S > 0 (as needed for non-singular rotating black holes with 
negative cosmological constant [14, 29]), to calculate the Hamiltonian mass we 
need to find the leading order behaviour of the metric and compare it to anti 
de Sitter. Eor this one needs first to transform the Boyer-Lindquist form of the 


K 


^ j dT^ + (^1 + -^ j + R^ © d<^^) (E.l) 
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metric to a new coordinate system defined by (see, e.g., [2]) 



= ^(^2^, + a2sin2(0)) , 

-Rcos( 0 ) = rcos( 0 ), 

^ A 

$ = 09 — a — t , 

^ 3S ’ 


Under (E.2) we have, with expansions for large R, 


(E.2) 


gTT 


qt^ 


9RR 


9Re 

gee 

g^^ 


^ gtt + 


2HaA 


3 

- 1 + 


9tip + 


a?h? 


9 


9<f«P 


72\/6; 


m 


R {a?K - a?k cos(20) + 6 )®/^ 


+ 0 (i? 


- 2 a 


^gtip + 2 gipip 


aA 

T' 

72 (\/6amsin^(0)) 


i? (a 2 A - a 2 Acos( 20 ) + 6 )®/^ 
3 9 


+ 0{R 


-2 a 


Ai?2 A2RA 


+ 


108\/6^ 


m 


A2i?5 (a 2 A - a 2 Acos( 20 ) + 6 )^/^ 
216y/6a^m sin( 20 ) 


+ 0 (i? 


-6a 


Ai?4 (a 2 A _ a 2 Acos( 20 ) + 6 )®/^ 

^2 lO8\/6a^msin^(20) 

7?3 (a 2 A _ a 2 Acos( 20 ) + 6 )^/^ 
sin2(0) 


+ 0 (i? 


-5a 


+ 0 (i? 


-4a 


+ - 


72\/6a^msin^(0) 


i? (a2A — a2Acos(20) + 6) 


5/2 


+ 0 (i? 


-2a 


(E.3) 

(E.4) 

(E.5) 

(E. 6 ) 

(E.7) 

(E. 8 ) 


To calculate the Hamiltonian mass of the Kerr-anti de Sitter metrics we can 
use the results of [21]: For this, let b = b^^dx^dx’^ denote the anti de Sitter 
metric in the coordinate system (4.29). Let e^, a G {0,1, 2, 3} be the following 
ON frame for b, 


1 + 




i?2 1 1 

dr, ej = ^/ 1 + -^dR , 62 = -de , eg = ■ (E.9) 


Define 

:= - b >^‘', (E.IO) 

and let e“^ denote the components of in the coframe dual to {ea}- If there 
exists e > 0 such that 

= 0(i?-3/2-^) ^ ^ ^ 
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then the Hamiltonian mass Mh of a hypersurface {t = r} equals [21, Equa¬ 
tion (5.22)] 


Mh = lim 




R-^oo 167r£2 Jg 


[ {RdRt 


eAA _ 2e^^)(fS 


= lim 
R-^oo 


/ ( - Rdne^j^ + 2eii)d S 

ru/ic J St,r 


(summation over A). 
One finds 


i?2 


"11 


= A + ^)i9RR - bRR) 

36 (\/6m) 


i?3 (A (a2A - a^A cos(20) + 6)^/^ 


+ 0{R 


-4a 


"12 


^22 ~ 


= OiR-^), 

{gee - bee) 


®33 ~ 


®oo ~ 


®03 ~ 


i?2 

108\/6o^m sin2(20) 

R^ (a^A — a2Acos(20) + 6)'^'^^ 

sm (0) 

72\/6a^msin^(0) 
i?3 (a2A - a2Acos(20) + 6)^/^ 

216 { V6m) 


R^ (a (a2A - a2A cos(20) + 6)^^^ 

216 {\/2am sin3(0)) 

R (^^/^ (a2A - a2Acos(20) + 6)®/^) 


+ 0{R-^), 

+ OiR-^), 

+ 0{R-^), 

+ 0(i?-2), 


(E.12) 


(E.13) 

(E.14) 

(E.15) 

(E.16) 

(E.17) 

(E.18) 


with vanishing remaining components. The fall-off requirements are therefore 
satisfied, and we obtain a mass integrand 

36 {s/Gm (a^(-A) -|- a^Acos(20) -|- 12)) 

-7-777\- ^(^{^ ) 

i?3 (a (a2A — a2Acos(20) -|- 6)^'^ j 

which integrates over 3“^ to 


Mh 


9m 

(a2A + 3)2 



(E.19) 


F Conformal Yano-Killing tensors 

In four-dimensional spacetimes admitting non-trivial conformal Yano-Killing 
(CYK) tensors Yag, or asymptotic CYK tensors, global invariants can be de- 
hned by integrating over two-dimensional submanifolds (cf., e.g.. 
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[31,32]and references therein). In Kerr-de Sitter spacetime [35] a solution of 
the CYK equations is given by 

Y := r sin 9d6 A [(r^ + a^) d(p — adt] + a cos 0dr A (dt — a sin^ 0d())) . (F.l) 

The Hodge dual of Y is also a CYK tensor: 

*Y = acos 0sin6*d0 A [(r^ + a^) d(/> — adt] + rdr A (asin^ 6<d(/> — dt) . (F.2) 

According to [3, p. 17], these two tensors form a basis of the set of solutions. 
As discussed extensively above, there is no unique choice of a preferred Killing 
vector X which can be used to define mass. Similarly, when we define the mass 
via CYK tensor, we have a freedom to multiply Y by any function of m and 
a, and there does not seem to be a preferred choice for asymptotically KdS 
metrics. 

The Weyl tensor for KdS depends on the cosmological constant A in a 
non-trivial way. Surprisingly enough, the CYK-Weyl contractions F{C,Y) := 
Y^’^Cf^uXn and F{C,*Y) := *Y^'^C^uXk do not depend on A, which results in 
the following closed two-forms: 


F{C,Y) = 

{[r^ — cos^ 6] sin 9d0 A [(r^ -|- a^)d(p — adt] 



+ 2ar cos 9dr A [a sin^ 9d9 — dt]} , 

(F.3) 

F{C,*Y) = 

1 2ar sin 9 cos 9d9 A [adt — (r^ + a^)d(p] 

P 



+ [r^ — cos^ 9]dr A (a sin^ 9d(p — dt)} , 

(F.4) 

where^ 

cos^ 9. Integration of (P.3) over a sphere r = 

const. 

and t = const, gives 


(F.5) 

while a similar integral for (F.4) vanishes. 
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